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We report results from a search for gravitational waves produced by perturbed intermediate mass 
black holes (IMBH) in data collected by LIGO and Virgo between 2005 and 2010. The search 
was sensitive to astrophysical sources that produced damped sinusoid gravitational wave signals, 
also known as ringdowns, with frequency 50 < /o/Hz < 2000 and decay timescale 0.0001 < r/s < 
0.1 characteristic of those produced in mergers of IMBH pairs. No significant gravitational wave 
candidate was detected. We report upper limits on the astrophysical coalescence rates of IMBHs 
with total binary mass 50 < M/M@ < 450 and component mass ratios of either 1:1 or 4:1. For 
systems with total mass 100 < M/ Mg < 150, we report a 90%-confidence upper limit on the rate 
of binary IMBH mergers with non-spinning and equal mass components of 6.9 x 10 -8 Mpc _3 yr _1 . 
We also report a rate upper limit for ringdown waveforms from perturbed IMBHs, radiating 1% of 
their mass as gravitational waves in the fundamental, I = m = 2, oscillation mode, that is nearly 
three orders of magnitude more stringent than previous results. 

PACS numbers: 95.85.Sz, 04.70.-s, 04.80.Nn, 07.05.Kf, 97.60.Lf, 97.80.-d 


I. INTRODUCTION 

Intermediate mass black hole (IMBH) binary systems 
represent a potential strong source of gravitational radia- 
tion accessible to ground-based interferometric detectors 
such as the Laser Interferometer Gravitational- Wave Ob- 
servatory (LIGO) [T] and Virgo [2]. Although yet to be 
discovered, binary systems with total masses in the range 
50 < M /Mq < 10 5 could form in dense star clusters such 
as globular clusters m- 

The coalescence of a compact binary system gener- 
ates a gravitational wave signal consisting of a low fre- 
quency inspiral phase when the compact objects are in 
orbit around each other, a merger phase marking the co- 
alescence of the objects and the peak gravitational wave 
emission, and a high frequency ringdown phase after the 
objects have formed a single perturbed black hole Ed- 
For low mass systems, most of the signal-to-noise ratio 
comes from the inspiral phase of the coalescence. Sev- 
eral searches for gravitational waves from the inspiral of 
low mass compact objects have been performed by LIGO 
and Virgo I®. However, since the merger frequency 


is inversely proportional to the mass of the system, it 
is shifted to lower frequencies for higher mass binaries. 
Searches for gravitational waves from the inspiral, merger 
and ringdown of binary black holes with total masses 
25 < M/ Mq < 100 have also been performed in LIGO- 
Virgo data 'mm2]. 

For an IMBH binary, typically only the merger and 
ringdown parts of the signal fall above the low frequency 
cutoff of 40 Hz for the initial LIGO and Virgo detectors. 
Thus it is sufficient to conduct a search solely for these 
particular phases of the gravitational wave signal [13l - 
ng. a binary black hole merger is expected to result in 
a single perturbed black hole, and black hole perturba- 
tion theory and numerical simulations provide us with 
a well-understood ringdown signal model, a superposi- 
tion of quasinormal modes that decay exponentially with 
time jTTTITSTTl . Indeed, any perturbed black hole, not just 
that produced by a compact merger (e.g., a black hole 
formed as the result of the core collapse of a very mas- 
sive star [MHMj), will emit ringdown gravitational waves 
described by its quasinormal modes. 

Since the gravitational waveform of perturbed black 
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holes lias a well-defined model, the method of matched 
filtering is used to search for ringdown signals. The first 
such search was carried out on data from the fourth LIGO 
science run (S4) which took place between February 22 
and March 24, 2005 [27] . Additionally, two burst searches 
with less-constrained waveform models looked for gravi- 
tational waves from mergers of IMBHs in data collected 
by LIGO and Virgo between 2005 and 2010 [251 125] . No 
events were observed in these searches. In this paper, we 
present the results of a matched filter ringdown search of 
data from LIGO’s fifth and sixth science runs and Virgo’s 
science runs 2 and 3. We compare the resulting rate up- 
per limits to the previous searches for gravitational waves 
from IMBHs. 

Sections m and m describe the expected ringdown 
sources and waveform. Section[ll]provides a brief descrip- 
tion of the detectors and their sensitivities during the 
data collection epochs. Section m describes the search, 
and results are presented in Section |TV} Upper limits are 
presented in Section [V] and discussed in Section VI 


A. Ringdown sources 

Observed black holes of known masses fall into two 
broad mass ranges. Stellar mass black holes have masses 
< 35 Mq [5U1452] although theoretical modeling of stel- 
lar evolution and population synthesis raises the possi- 
bility that significantly heavier stellar black holes could 
exist [33;, [31]. Supermassive black holes have masses 
> 10 5 M 0 [35J [35] and are thought to be cosmologi- 
cal in origin, possibly formed through galactic mergers 
leading to their growth through coalescences and accre- 
tion [3?l|33!. The large gap between the mass ranges of 
stellar and supermassive black holes is predicted to be 
populated by an elusive class of objects known as inter- 
mediate mass black holes (IMBHs) [3] 1394452] . Obser- 
vational evidence from ultra- or hyper-luminous X-ray 
sources and star cluster dynamics suggest a population 
of IMBHs with masses in the range 10 2 M 0 to 10 4 M 0 [3|. 
Ultra-luminous X-ray sources with angle-averaged fluxes 
many times that of a stellar mass black hole accreting at 
the Eddington limit (> 3 x 10 39 erg s _1 ) may be explained 
by black holes with masses larger than any known stellar 
mass black hole. The brightest known hyper-luminous 
X-ray source and the strongest IMBH candidate is the 
point-like X-ray source HLX-1. Its maximum X-ray lumi- 
nosity of 10 42 ergs _1 requires a black hole mass > a few 
10 3 M 0 [351 144], Other hyper-luminous X-ray sources 
include M82 X-l J35j, Cartwheel N10 [33], and CXO 
J122518.6 [47]. Furthermore, the excess of dark mass 
at the centers of globular clusters could be explained by 
~ 10 3 M 0 IMBHs formed from repeated mergers between 
other compact objects and/or stars [3544551 . However, 
both hyper-luminous X-ray sources and central globular 
cluster masses can be explained via phenomena that do 
not include IMBHs [53 1321 . Still, most observational ev- 
idence for globular cluster IMBHs using radio emissions 


can place upper bounds of < 10 3 M 0 [554455] , and do not 
rule out lower mass systems that are above the expected 
maximum mass of a normal stellar mass black hole [33] . 
Thus, the existence of IMBHs currently remains specu- 
lative. 

Numerical simulations suggest that IMBH binaries 
could form in collisional runaway scenarios in young 
dense star clusters. Initially, in young star clusters, 
IMBHs could form via the runaway collapse of very mas- 
sive stars [41] [591461] . After separate formation, two 
IMBHs could settle to the core of the cluster through 
dynamical friction and form a common binary via dy- 
namical interactions. The binary would tighten due to 
three-body encounters, finally merging quickly via grav- 
itational radiation [4, [52] [63] . 

From [53], we know that the astrophysical rate of 
IMBH binary coalescence in globular clusters (GC) 
should be no higher than 0.07 GC~ 4 Gyr _1 assuming that 
all globular clusters are sufficiently massive and have 
a sufficient binary fraction to form this type of binary 
once in their lifetime of 13.8 Gyr [5]. Also, globular clus- 
ters have a space density of roughly 3GCMpc -3 55] . 
This allows us to convert the astrophysical upper limit 
to 2 x 10 _1 ° Mpc _3 yr _1 . If we assume that only 10% of 
globular clusters meet these requirements, the rate would 
still be as high as one tenth this value [531 . 

Numerical simulations also suggest the possibility of 
forming intermediate mass ratio inspirals (IMRIs) (e.g., a 
coalescence of an IMBH with a compact stellar mass com- 
panion) in these same dense star clusters. This occurs 
through a combination of gravitational wave emission, 
binary exchange processes, and secular evolution of hi- 
erarchical triple systems [32 } [664469] . Ringdown searches 
in the advanced detector era could be important for de- 
tecting IMRIs, particularly if the inspiraling companion 
is a black hole with m > 10 M 0 or if the system is a 
compact object coalescing with a slowly-spinning IMBH 
with m > 350 M 0 [65] . 


B. Ringdown waveform 

A black hole can be perturbed in a variety of ways, 
e.g., by interaction with a companion, by accretion or 
infall of matter, or in its formation through asymmet- 
ric gravitational collapse. A perturbed Kerr black hole 
will emit gravitational waves, relaxing to a stable config- 
uration through radiation generated by a superposition 
of quasinormal modes of oscillation [364423] . The emit- 
ted gravitational waves are exponentially decaying sinu- 
soid signals characterized by a complex angular frequency 
Wfmn from which we can derive both the real frequency 
famn and the quality factor Q tmn : 

fgmn — ^(^mn)/2 7T , (1) 

Qtmn = Tt fimn/^S (S^Zran) 7 ( 2 ) 

where l = 2,3,..., and m = t are the spheroidal 

harmonic indices and n denotes the overtones of each 
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mode. Overtones with n > 0 are generally negligible 
in amplitude compared with the fundamental n = 0 
mode. Numerical simulations have demonstrated that 
the £ = m = 2 fundamental mode dominates the gravita- 
tional wave emission, particularly in the case of an equal 
mass compact object merger DS3- The ringdown search 
uses single-mode waveform templates. However, other 
modes can contribute significantly to the gravitational 
wave signal, particularly in cases where the binary’s mass 
ratio q = m>/m< ^ 1 where m> = max(mi,m 2 ) and 
m< = min (mi, m 2 ). Reference [7ll reports that single- 
mode templates can result in a loss > 10% in detected 
events over a significant mass range and also result in 
large errors in the estimated values of parameters (espe- 
cially the quality factor). A multimode ringdown search 
would perform better both in efficiency and parameter 
estimation m- Nevertheless, we show that the single- 
mode ringdown search will still provide good sensitivity 
to comparable mass binar y syst ems (see average sensitive 
distances given in Section VB). 


The response of an interferometric detector to a grav- 
itational wave is 


h(t) = F +{ e A,^)h+{i) + F x (6,<p,ip)h x {t) (3) 


square of the symmetric mass ratio, v = rn-i rn^ / (mi + 
m 2) 2 = g/(l + g) 2 , as e « 0.44V 2 [72U74] . Thus, for q = 1, 
e ~ 3% and, for q = 4, e ~ 1%. Gravitational waves from 
extreme mass ratio systems will not be detectable unless 
the system is sufficiently close (see Section VB| . A black 
hole binary with spinning components will radiate more 
energy if the spins are aligned with the orbital angular 
momentum and less if the spins are anti-aligned [73 [75]. 

The black hole mass M and dimensionless spin pa- 
rameter a can be determined numerically using fitting 
formulae to Kerr quasinormal mode frequency and qual- 
ity factor parameters tabulated in Table VIII of [TBj ■ For 
the {£, m, n ) = (2, 2, 0) mode, the fits are of the form: 


fo = 


1 


2t r GM L 


1.5251 - 1.1568(1 -a) 


0.1292 


Q = 0.7000 + 1.4187(1 -a) 


- 0.4990 


( 7 ) 

(8) 


These fitting functions allow us to relate a measurement 
of the frequency and quality factor from a match filter 
ringdown template to the mass and angular momentum 
of the final perturbed black hole. 

We can approximate the ringdown gravitational wave 
strain by 


where F + and F x are the antenna pattern functions that 
depend on the direction to the source as described by a 
polar angle 9, an azimuthal angle <j > , and a polarization 
angle ip. The plus and cross polarizations h+ and h x of a 
single-mode ( £,m,n ) = (2,2,0) ringdown waveform take 
the approximate form 

h + {t- 1 , 0) =y (1 + cos 2 i) e-*M‘-‘o)/Q 
cos [2tt/o (t - to) + (po ) , 


h 0 (t)=A e ge irfo{ ' t to)/Q cos[2Tif 0 (t-to) + tpo}, (9) 

for t > to where A c g = A/V e g and D e g is the effective 
distance to the source and ipo is the effective initial phase 
depending on the initial phase (po as well as on the signal 
polarization [see Eq. (1.7) and (1.9) in [77]]. Note that 
both (po and time of arrival at the detector to are set 
to zero for simplicity in the template waveform given in 
Section IIII Al 


h x {t; 1 , <p) = — (2 cos 1 ) e w 7o(* ( o )/Q 
r 

sin [2 tt/ 0 (t - t 0 ) + (p 0 \ , 


( 5 ) 


for t > to where fo = /220 and Q = Q 220 are the oscil- 
lation frequency and the quality factor of the (£, to, n) = 
(2, 2, 0) mode, r is the distance to the source, (po is the 
initial phase of the mode, and 1 is the inclination an- 
gle. The oscillation amplitude of the ( £,m,n ) = (2,2,0) 
mode, A, is given approximately by (see Appendix |A| 

_ 4 = gm /| 0 _ i, 2 ' - V2 , 

C z V 2 


where G is the gravitational constant, M is the black 
hole mass, c is the speed of light, e, known as the ring- 
down efficiency, is the fraction of the black hole’s mass 
radiated, a = cS/GM 2 where S is the black hole’s spin 
angular momentum, F(Q) = 1 + 1/(4Q 2 ) and q (a ) = 
[1.5251 - 1.1568(1 - a) 0 1292 ] [cf. Eq. (Frl, ®, and |A5|)]. 

The total ringdown efficiency of a blacK hole binary 
with non-spinning components is known to scale with the 


II. DATA SET 

The data analyzed spans multiple science runs for both 
the LIGO and Virgo detectors. We report results both for 
data collected between November 2005 and September 
2007 and between July 2009 and October 2010. 

The first time period covers LIGO’s fifth science run 
(S5). The LIGO site in Hanford, Washington hosted 
two collocated interferometers: a 4 km detector HI and 
a 2 km detector H2. The LIGO site in Livingston, LA 
hosted one 4km detector LI. Additionally, the Virgo 
3 km detector in Cascina, Italy operated from May 2007 
to September 2007 during its first science run (VSR1) 
which overlapped with the last few months of LIGO’s S5 
run. However, this search did not analyze VSR1 data. 
Thus, for the first time period, which we designate Pe- 
riod 1, we report results for the three-fold coincident 
search of the H1H2L1 detector network. We also report 
results for two-detector combinations of this network in- 
cluding H1L1 and H2L1. We chose to exclude H1H2 coin- 
cident events since accurately measuring the significance 
of gravitational wave candidates is complicated by this 
network’s correlated detector noise. 


The second time period covers LIGO’s sixth science 
run (S6) during which only the HI and LI LIGO detec- 
tors were operating. The Virgo detector conducted two 
science runs during this period: VSR2 which ran from 
July 2009 to January 2010, and VSR3 which ran from 
August 2010 to October 2010. For this second time pe- 
riod, which we designate Period 2, we report results for 
the coincident search of the H1L1V1 detector network. 
We also report results for all two-detector combinations 
within this network. 

LIGO’s S5 run marked the final data collection of 
the initial LIGO detector configuration during which de- 
sign sensitivity was achieved PP- Figure [l] (left) demon- 
strates the HI, H2, and LI detectors’ sensitivities to ring- 
down signals from spinning black holes with a = 0.9 and 
e = l/flj for typical Period 1 performance. This figure 
shows the horizon distance Uh divided by the square 
root of the ringdown efficiency e, scaled to a canonical 
value e = 1%, as a function of the final black hole mass. 
The horizon distance is the distance at which a given 
source with optimal location and orientation would pro- 
duce a SNR of 8 in a given detector; some details of its 
derivation for ringdowns are given in Appendix [B] Dips 
in the ringdown horizon distance correspond directly to 
features of the detectors’ noise spectral density curves. 
For instance, the strong dip in sensitivity at 360 M 0 is 
due to 60 Hz electric power noise. 

The S6 run, during the phase of the enhanced LIGO 
detector configuration, followed a series of upgrades to 
the initial detectors to improve sensitivity. These en- 
hancements included a higher power laser and a new DC 
readout system m- Similarly, the Virgo detector saw 
several improvements between its VSR1 and VSR2 runs 
including a more powerful laser, a thermal compensation 
system, and improved scattered light mitigation. Be- 
fore Virgo’s VSR3 run in early 2010, monolithic suspen- 
sions with fused-silica fibers were installed [79] . Figure [l] 
(right) demonstrates the HI, LI, and VI detectors’ sen- 
sitivities to ringdown signals from spinning black holes 
with a = 0.9 and e = 1% for typical Period 2 perfor- 
mance. 

Gravitational-wave strain data from each of the de- 
tectors are known to be both non-Gaussian and non- 
stationary. Non-Gaussianity is often manifested as noise 
transients, or glitches, in the strain data. Efforts are 
made to diagnose and remove glitches and stretches of 
elevated noise from the data set using environmental and 
instrumental monitors [80L182] . In this search, as in pre- 
vious searches of LIGO-Virgo data, we apply three levels 
of data quality vetoes [SB] SHI] (see Appendix A of [8] for 
more details). Data remaining after the first and second 
veto levels have been applied are searched for possible 
detection candidates (see Section IV I . Data remaining 
after all three veto levels have been applied are searched 


1 These values were chosen so that a direct comparison could be 
made with Fig. 2 in ED- 


TABLE I. Length of each network’s total analyzed time after 
the third level of vetoes has been applied and the playground 
data set has been removed. 


Analysis Tim^] (years) 

Network Period 1 Period 2 


H1L1 

0.09 

0.17 

H1V1 

- 

0.10 

H2L1 

0.07 

- 

L1V1 

- 

0.06 

H1H2L1 

0.63 

- 

H1L1V1 

- 

0.08 

Total 

0.79 

0.41 


a Excluding playground time. 


for detection candidates and are also used in constrain- 
ing the IMBH merger rate (see Section [V]). Table [I] gives 
the total analyzed time after all three veto levels are ap- 
plied and after the removal of the “playground” data set 
used for pipeline tuning as described in Section [TlI D The 
total analysis time for both Period 1 and Period 2 was 
1.2 years. 


III. RINGDOWN SEARCH 
A. Search Algorithm 


The ringdown search algorithm, first introduced in in 
13, is based on the optimal method for finding mod- 
eled signals buried in Gaussian noise, the matched fil- 
ter [ 55 ] . The data from multiple gravitational wave detec- 
tors are match filtered with single-mode ringdown tem- 
plates to test for the presence or absence of signals in the 
data. The output is a signal-to-noise ratio (SNR) time se- 
ries m from which local maxima above a pre-determined 
SNR threshold, called triggers, are retained for further 
analysis. Since the noise in the detector data is non- 
stationary and non-Gaussian, matched filtering alone is 
not enough to establish that a trigger is a gravitational 
wave signal. Since detector noise can often mimic the 
signal for which we are searching, additional tests are 
employed including detector coincidence and SNR con- 
sistency. We use a search pipeline similar to the ihope 
pipeline described in [ 55 ]. Here we summarize the main 
steps of the ringdown search pipeline. 

The data conditioning and segmentation is discussed 
in detail in [87] . Each segment of data is filtered using a 
bank of ringdown templates characterized by frequency 
/o and quality factor Q. Following EH, the template 
used in this search is 

h(t) = e cos(27 r/ 0 t) , 0 < t < f max (10) 

[cf. Eq. ([9])], with a length of 10 e-folding times, f max = 
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Final Black Hole Mass (M 0 ) 



FIG. 1. Ringdown horizon distances, Du, weighted by the square root of the ringdown efficiency, yA , as a function of final 
black hole mass for Period 1 (left) and Period 2 (right). Here we have set e = 1%. The dimensionless spin parameter is set to 
a = 0.9. For example, during Period 1, a ~ 200 Mg ringdown source with e = 1%, a = 0.9, and optimal location and orientation 
at a distance of ~ 530 Mpc would produce a signal-to-noise ratio of 8 in the HI detector. 


IOQ/tt/o^ 

The template bank is tiled in (/ 0 , Q)-space according 
to the analytic approximate metric computed assuming 
white detector noise as described in mi cm ins] so that no 
point in the parameter space has an overlap of less than 
97% with the nearest template)^] The template parame- 
ters cover a frequency band between 50 Hz and 2 kHz and 
quality factor in the physical range between 2 and 20. 
This corresponds roughly to masses in the range 10 M 0 
to 600 Mq, and spins in the range 0 to 0.99. A fixed bank 
of 616 templates was used for all detectors. 

Triggers with an SNR statistic above a predetermined 
threshold p* are retained for further analysis. For both 
Period 1 and Period 2, we set = p\ A = 5.5. For the 
least sensitive detector in each analysis period, we set 
lower thresholds: pi l2 = 4.0 and py 1 = 5.0. 


2 An arbitrary initial phase parameter (or equivalently, a quadratic 

sum of sine and cosine template outputs) could be implemented 
in the template waveform to reduce the fraction of power lost in 
the event of a pure sine wave signal. The problem is most acute 
for the detection of perturbed black holes with high frequency 
(/o 1000 Hz) and low dimensionless spin parameter (a < 0.6) 

where significant power is lost by using a cosine template ED- 
However, allowing an arbitrary phase would increase the noise 
level of the search. Furthermore, any ringdown signal would 
follow a preceding waveform and there is some arbitrariness in 
the division of one from the other. 

3 The template placement metric is derived using a sine template 
in ESI whereas a cosine template is used to filter the data. Op- 
timally, the metric derivation should account for initial phase 
dependence as derived in m ■ In the high Q limit, the sine and 
cosine metrics coincide. 


B. Coincidence and Vetoes 

Once triggers are found in a single detector, we apply a 
coincidence test, analogous to the one introduced in ED], 
to check for multi-detector parameter and arrival time 
consistency. In order to include information about time 
coincidence dt and template coincidence for d/o and dQ 
in a single coincidence test, we construct a 3D-metric [88] 
to calculate the distances ds 2 3 between two triggers in (/ 0 , 
Q, f)-space. The quantity (1 — ds 2 ) is a measure of nor- 
malized signal mismatch. To account for the finite travel 
time between non-collocated detectors, we minimize ds 2 
for each detector pair over a range of allowed time differ- 
ences. Only pairs of triggers for which ds 2 < ds 2 = 0.4 
are kept as coincident candidates. During times when 
three detectors are operating, triple coincident events are 
constructed from sets of three triggers if each trigger in 
the set passes the coincidence test with every other one. 
We also consider H1L1 coincidences in a H1H2L1 net- 
work. 

We also apply second and third level vetoes to seg- 
ments of poor data quality as described in [55]. Addi- 
tionally, for Period 1, we apply a number of amplitude 
consistency tests that exploit the coalignment of HI and 
H2 [55] . These tests allow us to apply cuts to reduce the 
background of false alarms. 

C. Ranking Events 

Finally, the pipeline ranks the coincidences and deter- 
mines significance. For this purpose, a detection statis- 
tic is designed to separate signal-like coincidences from 
noise-like coincidences. Given the large number of pa- 
rameters that describe multi-detector coincidences, we 
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employ a multivariate analysis using cuts on multiple pa- 
rameters to help in classifying coincidences as signals or 
false alarms: i.e. , a multivariate statistical classifier. The 
parameters provided to the classifier to aid in charac- 
terizing the multi-detector coincidences included single- 
detector SNRs and differences in time and template pa- 
rameters between detectors, recovered effective distances, 
composite SNR statistics^] the 3D-metric distance be- 
tween triggers and the metric coefficients as well as data 
quality information from the hierarchical veto method 
described in m- Additional details of these parameters 
will be described in a future paper. 

To perform the multivariate analysis, we use a machine 
learning algorithm known as a random forest of bagged 
decision trees [521 [53] . Similar techniques have been im- 
plemented for detecting gravitational-wave bursts [51] 
and cosmic strings [95]. The training of the classifier 
uses two sets of data: a collection of coincidences associ- 
ated with simulated signals and a collection of accidental 
coincidences that act as a proxy for the background. 

The simulated signal set is generated by adding 
software-generated gravitational waveforms to the data 
and running a separate search. The simulated waveforms, 
described in more detail in Section |IIID included both 
full coalescence IMBH merger signals and lone ringdown 
signals. 

The set of accidental coincidences is generated using 
the method of time-shifted data that takes advantage of 
the fact that a real signal will produce triggers in each 
detector that are coincident in time. The data streams of 
detectors are shifted in time with respect to one another 
by intervals longer than the light travel time between 
sites plus timing uncertainties, then a search for coinci- 
dences is performed. These time-shifted coincidences are 
then almost certainly due to noise. For Period 1, the LI 
data stream was shifted by multiples of 5 seconds relative 
to HI and H2 for a total of 100 time-shifted analyses; the 
HI and H2 data streams were not time-shifted relative 
to one another. For Period 2, the LI data stream was 
shifted by multiples of 5 seconds and the VI data stream 
was shifted by multiples of 10 seconds relative to HI for 
a total of 100 time-shifted analyses. 

The classifier assigns a likelihood ranking statistic £ 
to each coincidence. A high likelihood implies the co- 
incidence is signal-like; a low likelihood implies the co- 
incidence is noise-like. For each candidate, we need to 
be able to assign a significance to its likelihood ranking. 
This is done by mapping a false alarm rate (FAR) to a 
candidate’s rank in order to assess its significance. We 
count the number of false coincidences in the time-shifted 
searches, record their likelihood values, and determine 
the analysis time T k of all the time-shift searches for a 
particular experiment time (e.g., H1L1 coincidences in 
a H1L1V1 network, H1L1V1 coincidences in a H1L1V1 


4 Some details of the composite SNR statistics used for classifica- 
tion are given in [El- 


network, etc.). We perform this calculation separately 
for each type of coincidence in each of the different ex- 
periment times. Then, for each candidate in each exper- 
iment, we determine the FAR at its likelihood value £* 
with the expression: 


100 

E N k (C>C*) 

FAR = fc=1 ^ (11) 

where N k is the measured number of coincidences with 
£ > £* in the k th shifted analysis. We performed a 
total of 100 time-shifted analyses. Finally, we can rank 
candidates by their FARs across all types of experiment 
times into a combined ranking, known as combined FAR, 
for a single experiment time as described in detail in 1961 . 
The combined FAR is the final detection statistic that 
allows us to combine the candidate rankings from the 
various experiment types into a single list of candidates 
ordered from most significant to least significant. 


D. Tuning and simulations 

The analysis was tuned using the set of false alarm 
coincidences obtained from time-shifted searches, a set 
of simulated signals (“injections”) added to the detec- 
tors’ data streams in a separate stage of data analysis, 
and a small chunk of the actual search data, approx- 
imately 10%, designated “playground”, that was later 
excluded from the analysis to preserve blindness. The 
goal of tuning the analysis is to maximize the sensitiv- 
ity of the search while minimizing the false alarm rate. 
For this, we injected a set of ringdown-only waveforms 
with e = 1% into the data set. The waveforms were 
determined by Eq. 0 , 0 , and 0 with sky location 
and source orientation sampled from an isotropic distri- 
bution. Several sets of ringdown waveforms were injected 
with a uniform distribution in /o and Q to cover the pa- 
rameter range of the ringdown template bank. Also, in 
order to cover the broad mass and spin range accessi- 
ble to the ringdown search when signals have e = 1%, 
several sets of ringdown waveforms were injected with a 
uniform distribution in M and d: 50 < M/ M Q < 900 and 
0.0 < d < 0.99. Additionally, we also injected a set of full 
coalescence waveforms with isotropically-distributed sky 
location and source orientation parameters into the data. 
These full coalescence waveforms included the recently- 
implemented non-spinning EOBNRv2 family m and the 
spinning PhenomB family [95]. The EOBNRv2 injections 
were distributed uniformly in total mass 50 < Mj Mg < 
450 and in mass ratio 1 < q < 10. The PhenomB injec- 
tions were given the same mass distribution and a uni- 
form dimensionless spin parameter 0.0 < 01,2 < 0.85 
where dp 2 = cSi^/Gml 2 for the spin angular momen- 
tum S and the mass m of the two binary components. 
For a discussion of the injection sets used in computing 
rate upper limits, see Section [V] 
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IV. SEARCH RESULTS 

The search yielded no significant gravitational wave 
candidates, as all events were consistent, within 1 sigma, 
with the background from accidental coincidences. Fig- 
ure [2] shows the cumulative distributions of coincident 
events found as a function of inverse combined false alarm 
rate after all vetoes up to the third level are applied. 
These plots combine results from both triple and double 
coincident searches over the total analysis time of Pe- 
riod 1 and Period 2. 

The most significant event was found in triple coin- 
cidence during Period 1 in HI, H2, and LI. After the 
first and second level vetoes were applied, it was found 
with a combined FAR = 2.07 yr -1 and, after the third 
level vetoes were additionally applied, with a combined 
FAR = 0.45 yr -1 . Thus we expect an accidental coin- 
cidence to be found by the search with this significance 
~ once per two years of analysis. Since the total anal- 
ysis time was 1.2 years, the event is consistent, within 1 
sigma, with the accidental coincidence rate. In both HI 
and H2, a trigger was found barely above threshold with 
matched filter SNRs of 5.5 and 4.4, respectively. How- 
ever, the candidate was found as a very loud trigger in 
LI with a matched filter SNR of 48.9. Performing a co- 
herent Bayesian parameter estimation follow-up [99] on 
these triggers, we found that a coherent analysis favored 
a solution for the binary’s sky location and orientation 
that yield a very strong signal in LI, but virtually no re- 
sponse in HI and H2 detectors. While it is theoretically 
possible that very particular location and orientation pa- 
rameters could produce such a signal, an excursion from 
stationary, Gaussian noise (a glitch) in LI is more likely. 


V. RATE LIMITS 

In this section, we compute the 90%-confidence upper 
limits on IMBH coalescence rates and IMBH black hole 
ringdown rates. The former will allow us to make an 
astrophysical statement as well as to compare the sen- 
sitivity of the ringdown search to various other searches 
that have made statements in this mass regime, includ- 
ing [nmaEsiEg. 

We used a procedure similar to that discussed in El 
Q1 for the upper limit calculation based on the loudest 
event statistic rmni mm . In order to capture the vari- 
ability of the detector noise and sensitivity, we analyzed 
the data in periods of ~ 1 to 2 months. In each of these 
analysis times, we estimate the volume to which the ring- 
down search is sensitive by injecting many simulated sig- 
nals into the data and performing an analysis to recover 
them. In Section |V B[ we describe the distribution of 
EOBNRv2 waveforms used to model the source popula- 
tion of IMBH binaries. Our sensitivity to these signals 
depends on total mass, mass ratio, source distance, and 
sky location as well as other parameters such as compo- 
nent spins. We explore the changing sensitivity of the 



FIG. 2. Cumulative distributions of coincident events found 
as a function of inverse combined false alarm rate after all 
vetoes up to the third level are applied. The figures combines 
results from both triple and double coincident searches over 
the total analysis time of Period 1 and Period 2. Grey con- 
tours mark the la through 5a region of the expected back- 
ground from accidental coincidences. No search candidates 
stand out from the background. 


ringdown search to these binaries over a range of total 
masses for both equal mass and 4:1 mass ratio systems. 
Other distance and orientation parameters are randomly 
sampled. Due to the significant variation of the search 
sensitivity over the large mass and mass ratio parameter 
space that we explore in Section VB[ we have chosen to 
include only systems with non-spinning components in 
this study. In Section |V C[ we describe the distribution 
of ringdown waveforms used to model the population of 
perturbed black holes first explored in [27] . 

For each of these injection sets, we compute the sen- 
sitive volume for a given mass range and mass ratio by 
integrating the efficiency of the search over distance: 


V eS = 47 t J r/(r)r 2 dr 


( 12 ) 


where the efficiency 77 (r) is calculated as the number of in- 
jections found with a lower combined FAR than the most 
significant coincident event in each analysis time for the 
search divided by the total number of injections made 
at a given distance. As described in [m H 21 nooi noT| . 
we estimate the likelihood parameter A of the loudest 
event being a signal versus being caused by an acciden- 
tal coincidence for each type of coincident network time 
and each mass and mass ratio bin. For each analysis 
time (excluding playground time), effective volume from 
Eq. (12), and estimated A, we marginalize over statis- 
tical uncertainties given in Section |V A| and construct a 
marginalized likelihood as a function of the astrophysi- 
cal rate in units of mergers per Mpc 3 per year for our 
EOBNRv2 injection sets and in units of ringdowns per 
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Mpc 3 per year for our ringdown injections. In order to 
obtain a combined posterior probability distribution for 
the rate over all the analysis times, we multiply a prior 
on the rate by the product of the marginalized likelihood 
functions to obtain a posterior probability and integrate 
to 90% to obtain the 90%-confidence upper limit on the 
rates. For our combined Period 1 result, we assumed a 
uniform prior on the rate. However, for the main Pe- 
riod 2 result, we were able to use the Period 1 posteriors 
over coalescence or ringdown rate as priors for the upper 
limit calculation. 


A. Sources of uncertainty 

We must account for several sources of random and 
systematic error when computing rate upper limits. Un- 
certainties on the sensitive volume as well as incomplete 
knowledge of waveforms and source populations form the 
largest contributors. As described in earlier search pa- 
pers )10H12| . we marginalize over random uncertainty 
(i.e. calibration and statistical Monte Carlo uncertain- 
ties) for each analysis time. The 90%-confidence upper 
limits based on the marginalized posterior distributions 
are the main results of this search. 

The calibration of the data is a source of both random 
and systematic error. Reference [102] reports uncertain- 
ties on the magnitude of the response function for each 
detector in Period I. We find an overall distance uncer- 
tainty of 8%. Thus, the random uncertainty on the visible 
volume for Period 1 is approximately 8% cubed, or 24%. 
For Period 2, references 11031 and |l04j report uncertainty 
on h(t) for LIGO and Virgo detectors. Additionally, an 
uncertainty on the scaling of h(t) was reported in 1 103 1 
and should be treated as a systematic error similar to 
the systematic waveform uncertainties discussed below 
that could over- or under-bias the amplitude of a signal. 
However, the uncertainty on the scaling of h(t) also has 
an associated random error that we fold into the random 
uncertainty calculation for Period 2. We find an overall 
distance uncertainty of 14% corresponding to a 42% un- 
certainty on the visible volume for Period 2. See 1 105] 
for a detailed explanation of how the uncertainties were 
propagated. 

In addition to the systematic error associated with the 
overall scaling of h(t) that could lead to amplitude bias 
as mentioned above, there is a larger source of systematic 
error due to differences between the injected model wave- 
forms and the true waveform. For EOBNRv2 waveforms 
below ~ 250 M 0 , comparisons with numerical models 
indicate that uncertainties in these waveforms result in 
< 10% systematic uncertainty in the SNR, corresponding 
to a < 30% uncertainty in sensitive volume. For higher 
masses, the systematic uncertainty in the SNR could be 
as high as 25%. Due to our incomplete knowledge of the 
true waveform and its changing uncertainty over the mass 
range we have explored, no systematic errors associated 
with imperfect waveform modeling were applied to the 


rate upper limits reported in this paper. Systematic er- 
rors were also not applied to previous searches mm us- 
ing full coalescence waveforms up to 100 M 0 and thus we 
can compare the upper limits directly with those results. 
A previous weakly modeled burst search [2B; used wave- 
form errors of ~ 15%. Thus, in order to compare with 
these results, the upper limits reported here should be 
rescaled as described below. Regarding ringdown wave- 
forms, due to our lack of knowledge about the population 
of black holes producing the waveforms and the wave- 
forms themselves, we again assign no systematic error to 
rate upper limits computed with ringdown waveforms. 

In general, we can rescale our rate upper limits by any 
systematic uncertainty by applying the scaling factor (1 — 
<t ) -3 where a is the systematic uncertainty. Thus, we 
can apply a conservative systematic uncertainty of 15% 
by rescaling our rate upper limit upward by a factor of 
1.63. 

The statistical error originating from the finite number 
of Monte Carlo injections that we have performed is the 
final source of error for which we must account. These 
errors on the efficiency at a given distance are found to 
range between 1.7% and 6.2% and were marginalized over 
using the method described in jlOO , 101] . 


B. Rate limits from full coalescence injections 

In order to evaluate the sensitivity of the ringdown 
search to waveforms from binary IMBH coalescing sys- 
tems with non-spinning components, we used a set of in- 
jections from the EOBNRv2 waveform family described 
in Section [ill D| Due to the variation in ringdown search 
sensitivity over different mass ratios, we chose to compute 
IMBH coalescence rate upper limits separately for <7=1 
and q = 4. The injection sets were distributed uniformly 
over a total binary mass range from 50 < M/ M 0 < 450 
and upper limits were computed in mass bins of width 
50 Mq. The final black hole spins of these injections can 
be determined from the mass ratios and zero initial com- 
ponent spins [I06j . For <7 = 1, we find a = 0.69, and for 
q = 4, we find a = 0.47. 

The average sensitive distances of the ringdown search 
to IMBH binaries described by EOBNRv2 signal wave- 
forms for both q = 1 and q = 4 are shown in Fig. [3] for 
Period 1 and Period 2. The most sensitive mass bin in 
both cases is 100 < Mj M 0 < 150 corresponding roughly 
to 110 < /o/Hz < 170 near the peak sensitivity of the 
LIGO detectors. For q = 1, the average sensitive distance 
of the 100 < M/Mq < 150 mass bin was 240 Mpc. For 
q = 4, the average sensitive distance for this mass bin 
decreases by more than a factor of two to 110 Mpc. As 
discussed in Section [lB| the reduced ringdown efficiency 
for q = 4 binary systems leads to lower amplitude wave- 
forms and hence, to lower average sensitive distances. 
Additionally, the lower final black hole spin for q = 4 
binary systems acts to decrease the average sensitive dis- 
tance relative to q = 1 binary systems for which the final 
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FIG. 3. Average sensitive distances of the ringdown search to 
binary systems described by EOBNRv2 signal waveforms over 
a range of total binary masses for Period 1 [q — 1 (yellow), 
q = 4 (green)] and Period 2 [q = 1 (cyan), q = 4 (blue)]. These 
distances are equivalent to appropriate averages over each of 
the detector networks shown for Period 1 and Period 2 in Ta- 
ble |T] weighted by the percentage of time analyzed for each 
network. Thus, while in general the H1L1V1 and H1L1 net- 
works during Period 2 were more sensitive than the H1H2L1 
and H1L1 networks during Period 1, the consistently smaller 
average sensitive distances for Period 2 reflect the large duty 
cycle of its least sensitive detector networks compared to Pe- 
riod 1 . 


spin is larger. The sensitive distance of higher mass bins 
drops off significantly due to the steeply rising seismic 
noise in the detector at low frequencies. This affect is ac- 
centuated for q = 4 systems relative to q = 1 systems at 
a fixed mass because a smaller final spin leads to a lower 
frequency ringdown. The sensitive distance of mass bin 
400 < AI/Mq < 450 is over an order of magnitude less 
than the sensitive distance of our most sensitive mass 
bins for both q = 1 and q = 4 cases. 

Figure [d] shows the 90%-confidence upper limits on 
non-spinning IMBH coalescence rates for a number of 
mass bins. We find an upper limit of 0.069xl0~ 6 Mpc~ 3 
yr -1 on the coalescence rate of equal mass IMBH bi- 
naries with non-spinning components and total masses 
100 < M/ M 0 < 150. From the discussion of astrophys- 
ical rates of IMBH mergers in Section EH we see that 
this rate upper limit is still several orders of magnitude 
away from constraining the astrophysical rate from GCs. 

Previous searches for weakly-modeled burst signals 
found no plausible events [251 [25] . The most recent search 
reports a rate upper limit for non-spinning IMBH coales- 
cences of 0.12 x 10 -6 Mpc _3 yr _1 at the 90%-confidence 
level for the mass bin centered on m\ = m 2 = 88 Mg [29] . 
A direct comparison of our q = 1 upper limits shown in 
Fig.[4]to this burst search result should be made with care 
due to the following differences between the two anal- 
yses: statistical approaches leading to different search 
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FIG. 4. Upper limits (90% confidence) on IMBH coalescence 
rate in units of Mpc~ 3 yr _1 as a function of total binary 
masses, evaluated using EOBNRv2 waveforms with 5=1 
(slate grey) and q = 4 (grey). In both cases, upper limits 
computed using Period 2 with Period 1 as a prior are shown 
in a darker shade. Overlaid in a lighter shade are upper limits 
computed using only Period 1 data with a uniform prior on 
rate. 


thresholds, treatment of uncertainties, analyzed detector 
networks, and mass and distance binnings. Additionally, 
while the ringdown search employed the Bayesian formu- 
lation flOO . IDT] for calculating the rate upper limit, the 
burst search used a frequentist method. Nevertheless, al- 
though the impact of the reported differences is hard to 
quantify, the upper limits determined by the two analy- 
ses can be considered consistent with each other. A more 
robust comparison of the sensitivity of the burst searches 
and an earlier version of the ringdown search without a 
multivariate classifier will be presented in a future pa- 
per IWl . 

Additionally, we can make a comparison with the 
upper limits reported from the matched filter search 
for gravitational waves from the inspiral, merger, and 
ringdown of non-spinning binary black holes with to- 
tal masses 25 < M/Mg < 100 jT2]. This search con- 
sidered similar uncertainties and similar analyzed net- 
works to those used by the ringdown search so a re- 
sult comparison is fairly straight-forward. From Table I 
of [H], we find that for systems with q = 1, the rate 
upper limits for masses 46 Mg to 100 Mg vary in the 
range 0.33x 10" 6 Mpc" 3 yr- 1 to 0.070x 10" 6 Mpc- 3 yr" 1 . 
From Fig. [4] we find a rate upper limit for mass bin 
50 < M/Mg < 100 of 0.16 x 10 -6 Mpc _3 yr _1 , a value 
consistent with the BBH rate upper limit range for these 
masses and mass ratio. 

Note that we can rescale our rate upper limits by 
a 15% systematic uncertainty by applying the scal- 
ing factor of 1.63 as described in Section |V A| From 
Fig- 0 we find a rescaled rate upper limit of 0.11 x 
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10 6 Mpc 3 yr 1 for mass bin 100 < M/ Mg < 150 and 
0.15x 10 -6 Mpc _3 yr _1 for mass bin 150 < Mj M 0 < 200. 


C. Rate limits from ringdown injections 


in Eq. Q and ([8|, differences in the method used in the 
volume integral in Eq. (12 1, and differences in marginal- 
ization over errors. A careful study of the improvement 
due to the use of a machine-learning algorithm will be 
presented in a future paper. 


In order to compare with E3, we determined a 90%- 
confidence upper limit of 4 x 10 -8 Mpc~ 3 yr _1 on rates 
of pure ringdowns from perturbed black holes with uni- 
formly distributed masses 85 < M/M 0 < 146, uni- 
formly distributed spins 0 < a < 0.99, and a fixed ring- 
down efficiency of e = 1%. We expect ringdown sig- 
nals from IMBH mergers to emit near this efficiency in 
the {JL = m = 2) fundamental mode if the mass ratio is 
near unity. However, for other sources of perturbed black 
holes, such as a hypermassive star collapse directly to a 
perturbed IMBH, we expect e -C 1%. Thus, the rate up- 
per limit reported in this section will not be applicable 
to such sources. 

Reference m placed a 90% confidence upper limit 
on the rate of ringdowns from black holes with frequen- 
cies distributed uniformly in logio (/o) in the range 70 < 
/o/Hz < 140 and uniformly in quality factor 2 < Q < 20 
of 3.2 x 10 -5 Mpc _3 yr~ 1 . Thus, a rough comparison in- 
dicates an improvement of nearly three orders of magni- 
tude. A significant portion of this improvement results 
from a huge increase in the analysis time. Due to the 
high false alarm rate in double coincident analysis time, 
an upper limit was set in m using only triple coinci- 
dent time, a total of 0.0375 years. We analyzed both 
triple and double coincident time in both Period 1 and 
Period 2, a total of 1.2 years. Such an increase in analysis 
time results in a factor of ~ 32 improvement in the upper 
limit. Additionally, a significant improvement in detec- 
tor sensitivity due to detector upgrades between science 
runs contributed to a better upper limit. Furthermore, 
since only triple coincident time was analyzed in [271 , the 
sensitivity was limited by the least sensitive detector, H2, 
which was shown to have a horizon distance of ~ 130 Mpc 
at 250 M 0 as shown in Fig. 2 in m- However, since we 
analyzed both triple and double coincident triggers, the 
limiting detector was typically the LI detector. We can 
compare the H2 horizon distance in Fig. 2 in m to the 
LI horizon distance in Fig. [l] at 250 M 0 to see that the 
horizon distance of the limiting detector improved by a 
factor of ~ 3 for a = 0.9. Since the upper limit scales 
with volume, a factor of ~ 3 in distance results in a fac- 
tor of ~ 27 in the upper limit. However, we expect this 
factor of improvement to decrease for the lower masses 
on which the ringdown upper limit was set. 

Thus, from the improvements both in analysis time 
and detector sensitivity, we find already roughly three 
orders of magnitude improvement. However, several 
caveats would apply to a direct comparison: different in- 
jection distributions in (M, a)-space, the improvements 
from pipeline enhancements such as the implementation 
of a machine-learning algorithm, differences in the fit- 
ting functions for final black hole mass and spin defined 


VI. SUMMARY 

This paper presents the results of the search for ring- 
down gravitational waves in data collected by LIGO and 
Virgo between 2005 and 2010. No significant gravita- 
tional wave candidate was identified. We place upper 
limits on the merger rates of non-spinning IMBH bina- 
ries as well as on the rates of ringdowns from perturbed 
black holes. 

We conducted a detailed study of the pipeline’s sen- 
sitivity to full coalescence IMBH merger signals using 
non-spinning EOBNRv2 waveforms. For simplicity, we 
focused our studies on only two mass ratios: q = 1 and 
q = 4. The average sensitive distances in our most sen- 
sitive total mass bin, 100 < M/ M 0 < 150, indicate that 
the ringdown search is sensitive to an equal mass system 
at twice the distance of a 4:1 mass ratio system. The most 
efficiently detected mass bin gives an upper limit on the 
rate of non-spinning, equal mass IMBH mergers with to- 
tal masses 100 < M /M 0 < 150 of 6.9 x 10 -8 Mpc _3 yr _1 . 
This does not account for any uncertainty in the wave- 
form, which could be as high as 10% for the mass bin. 
Our upper limits for ringdown waveforms from perturbed 
IMBHs with masses 85 < M/ M 0 < 146 and spins 
0 < a < 0.99 show an improvement of nearly three or- 
ders of magnitude over the previous result reported [27] , 
which we can attribute to improved detector sensitivity, 
increased livetime, and pipeline enhancements. 

While our rate upper limits are still two to three or- 
ders of magnitude away from constraining the astrophys- 
ical IMBH merger rate from globular clusters, we note 
that we will soon approach this optimistic rate with the 
improved sensitivity of Advanced LIGO and Virgo detec- 
tors expected to begin operation in 2015. With the im- 
proved low frequency performance of the advanced detec- 
tors, we will have sensitivity to gravitational waves from 
perturbed intermediate mass black holes with masses up 
to ~ 1000 to 2000 M 0 . At peak sensitivity, the Advanced 
LIGO ringdown horizon distance for black holes with 
e = 1% will approach cosmological distances. 
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Integrating this flux over a sphere with area element 
d A = r 2 d(cos i)d0, we find that the trigonometric func- 
tions simplify greatly, leaving only the exponential time 
dependence over which to integrate 


E = 


16ttG . 
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Finally, we note that the energy radiated as gravitational 
waves during the ringdown phase is E = eMc 2 where e 
is the ringdown efficiency discussed in Sec. m Thus, the 
amplitude can be found by solving Eq. ( A4 ) for A , 
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Appendix A: Ringdown Amplitude 


Appendix B: Ringdown Horizon Distance 


The amount of energy d E carried by gravitational ra- 
diation crossing an area dA orthogonal to its propagation 
direction in a time df is given by the energy flux equation, 


d E 
dAdt 


167tG 



(Al) 


where and h x are given by the generalized forms of 
Eq. © and © for an arbitrary location on a 2-sphere 
with to = 2 and time of arrival to set to zero. Taking the 
time derivative and squaring the plus and cross polariza- 
tions, we find 


h\ = (^)“(l + cos 2 0 2 e- 2 ^ 

(27t/ 0 ) 2 sin 2 (27r/ 0 t + 2 0) 

+ (“^) c °s 2 (27r/ 0 t + 2(f>) 

+(4tt/o) sin ( 27r /oi + 2</>) cos(2n f 0 t + 20) , 

(A2) 


The ringdown horizon distance, similar to the inspiral 
horizon distance, is a useful measure of the sensitivity 
of the detectors to ringdown gravitational waves from a 
particular type of black hole. It is equal to the distance 
at which an optimally oriented and located IMBH merger 
would produce an SNR of 8 in the detector. The horizon 
distance is derived from the representative strain noise 
power spectral density of a detector and the h ISS , or root 
sum squared of the strain, for a signal with optimal ori- 
entation at 1 Mpc. The definition of h rss comes from the 
need to measure the amplitude of a gravitational wave 
without reference to a particular detector. In general, it 
is 

POO 

h 2 ISS = / (h 2 + (t) + h 2 x (t)) d t, (Bl) 

Jo 

where h + and h x are given in Eq. Q and ([5]) for 
the single-mode (/, to, n) = (2, 2, 0) ringdown waveform. 
Here, under the assumption of optimal orientation, we 
set l = 0. We find that the h ISS takes the form 




e -2 nfot/Q 


(2tt/ 0 ) 2 cos 2 (2t xf 0 t + 20) 


where A is derived in Eq. (A5|. If h(f) represents the 
Fourier transform of the expected signal, then the average 
SNR this signal would attain in a detector with spectral 
density S n (f) is given by 


+ 



2 

sin 2 (2irfot + 20) 


-(4tt/ 0 ) 



cos(27r/ 0 f + 20) sin(27r/ 0 t + 20) . 

(A3) 



Typically, the horizon distance is found by setting (p) = 8 
and solving for the distance r which parameterizes the 


16 


waveform h. We can use the fact that the single-mode 
ringdown signal is quasi- monochromatic and S n (f) as- 
sumes approximately one value for each f 0 so it can be 
treated as a constant: 


(P) 


/ 4 f°° 


S n (fo)J 0 
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df. 


(B4) 


Using Parseval’s theorem, we can write Eq. ( B4 ) as 



Also, since optimally oriented and located sources imply 
maximization over all the angles 0, 0, and if) in F+ and 
E x , then F + = 1 and F x = 0. This then gives us the 
result that h(t) = h + {t) (which is defined for t > 0) so 


Eq. (B5) becomes 


(P) = 


S„(fo) 


h?±_(t)dt 


2 1 + 2Q 2 

:h 2 


Sn(fo) rss l + 4Q 2 


Sn(fo) 


ft L( 1M P c ) 


lMpcV 1 + 2Q 2 


(B6) 


) 1 + 4 Q 2 


where /i 2 ss (lMpc) is Eq. (B2) evalu ated at a distance of 
1 Mpc. Then, we simply solve Eq. (B6) for the horizon 
distance, 


r = 


1 Mpc 

(P) 


■I 


2 

S n (f 0 ) 


ft L( 1M P c ) 


1 + 2Q 2 

1 + 4Q 2 ' 


(B7) 


We then set (p) = 8 to define the ringdown horizon dis- 
tance used in Fig. [I] 
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